Zeeman effect in superconducting two-leg ladders: irrational magnetization plateaus 

and exceeding the Pauli limit. 
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The effect of a parallel magnetic field on superconducting two-leg ladders is investigated numeri- 
cally. The magnetization curve displays an irrational plateau at a magnetization equal to the hole 
density. Remarkably, its stability is fundamentally connected to the existence of a well-known mag- 
netic resonant mode. Once the zero-field spin gap is suppressed by the field, pairs acquire a finite 
momentum characteristic of a Fulde-Ferrell-Larkin-Ovchinnikov phase. In addition, S z = triplet 
superconducting correlations coexist with singlet ones above the irrational plateau. This provides a 
simple mechanism in which the Pauli limit is exceeded as suggested by recent experiments. 

PACS numbers: 71.10.Pm, 75.60.-d, 74.81.-g, 75.40.Mg 



Ladder systems, which consists of two or more strongly 
coupled chains, have been widely studied because of 
their exotic properties In particular, two-leg lad- 
ders are Mott insulators at half filling, with a spin gap, 
and become superconducting under doping (and pres- 
sure) 2]. This doped spin- liquid is believed to be a 
paradigm for the Resonating Valence Bond mechanism 
of superconductivity proposed by Anderson Such 
a low-dimensional superconductor provides an appeal- 
ing situation to study the role of a parallel magnetic 
field (Zeeman effect). In a singlet superconductor, the 
critical magnetic field which destroys pairing by split- 
ting up the up and down spin Fermi surfaces is called 
Pauli limit. The possibility of exceeding this theoreti- 
cal limit and/or stabilizing inhomogcneous superconduc- 
tivity predicted by Fulde, Ferrell, Larkin and Ovchin- 
nikov (FFLO) 0," have attracted strong interest in low- 
dimensional superconductors [jj. In particular, recent 
experiments suggest that Pauli limit is exceeded in the 
ladder material Sri4- :c Ca :!; Cu2404i (SCCO) 0,0. Fur- 
thermore, under magnetic field, insulating ladders can 
have plateaus for rational values of the magnetization, de- 
pending on the number of legs and on the interactions || . 
Doped low-dimensional strongly correlated systems such 
as doped spin chains 0, ^| have also proved to have in- 
teresting magnetization curves with plateaus at irrational 
values controlled by hole doping. Recently, Cabra et 
al. have predicted such magnetization plateaus in 
doped Hubbard ladders by means of Abelian bosoniza- 
tion and a strong-coupling expansion. However, numeri- 
cal evidence supporting such a scenario is lacking. 

In this letter, we investigate doped two-leg ladders with 
a Zeeman coupling u sing the density matrix renormaliza- 
tion group (DMRG) [12111.1]. We show that, for isotropic 
couplings givi ng rise to a superconducting Luther-Emery 
phase [bj, Il5| , a non-trivial plateau occurs for the finite 
magnetization value predicted in Yamanaka-Oshikawa- 
Affleck (YOA) theorem 0. We give a simple physical 
argument to explain this behavior which is similar to the 



formation of the resonant m agne tic mode encountered in 
the study of spin dynamics [lj|. We finally discuss the 
evolution of the ground state properties under magnetic 
field. In particular, we focus on the persistence of su- 
perconductivity for magnetic fields larger than the Pauli 
limit. This persistence is associated with the emergence 
of Cooper pairs with a finite total momentum, a typical 
feature of FFLO phases. We also observed the coexis- 
tence of singlet and triplet pairing. 

We describe the doped two-leg ladder using the t — J 
model with the same couplings along the legs and be- 
tween the legs. This model is believed to account for the 
unconventional superconductivity in lightly doped lad- 
ders, like the one observed in SCCO. We consider a static 
and uniform magnetic field H which couples to the spin 
degree of freedom via a Zeeman term. This corresponds 
to a magnetic field oriented parallel to the ladder plane to 
avoid orbital effect. Hence, we can write the Hamiltonian 

n l X 'PG[cl s c J , s + c[ s c hS \v G 

+J ^ I s * ' s i - -^iUj] - giL B y^HS, 

(id) ' l 

in which Vg is the Gutzwiller projector, c| s , Sj and rij 
are respectively electron creation, spin and density oper- 
ators at site i and s is the spin index. In what follows, 
the coupling constant g/is is absorbed in the definition 
of H. We computed energies at fixed magnetization us- 
ing the single center site method recently proposed by 
one of us which enabled us to have a discarded weight of 
order 10 -7 — 10 -8 with 1600 kept states and a noise level 
of 10~ 6 0. Energies are denoted by E(n hl S z ) with 
rih the number of holes and S z the total spin along the 
field. Open boundary conditions are used, L stands for 
the ladder length and S = nh/N s i te s is the hole doping 
(-^Vsitcs = 2L). All data were computed with J/t = 0.5 
for which the model is known to have dominant super- 
conducting correlations at low doping [bl. Hif . 
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FIG. 1: (Color online) Magnetization curves at fixed hole den- 
sity 5 — 0.0625 for different system sizes. Three plateaus are 
visible for m = 0, 5 and 1 — 5. H c is the superconducting criti- 
cal field (see text) and results in a discontinuity of the slope in 
the magnetization curve. The undoped case is shown for com- 
parison. Insert : Finite-size scaling of the two critical fields 
corresponding to the boundaries of the m = 5 plateau for var- 
ious densities 6 = 0.0625(0), 0.125(D), 0.1875(0), 0.25(A). 



Yamanaka-Oshikawa- Affleck theorem - We first recall 
the topological constraints that govern the opening of 
gaps in the excitation spectrum of one-dimensional like 
systems In the case of doped two- leg ladders, YOA 
relation takes the simple form 

l-<5±meZ, 

in which m = 2S Z /N s a es is the total magnetization nor- 
malized so that it equals 1 — 8 at saturation. It is a 
commensurability condition for either spin | or j but it 
introduces irrationality through 6. For 5 small enough, 
the relation can be rewritten 



= 8 



= 1-6. 



The theorem states that: (i) if the magnetization-doping 
relation is not satisfied, the low-energy spectrum can be 
either a continuum or gapped, but this latest case is pos- 
sible only if 1 — 8 ± m = p/q G Q and is associated 
with a spontaneous breaking of the translational symme- 
try so that the ground state is g-fold degenerate, (m) if 
the relation is satisfied, a gap can open in the excitation 
spectrum. Note that for doped systems with charge and 
spin degrees of freedom, one sector can be gapped while 
the other remains gapless. This actually happens for su- 
perconducting two-leg ladders: the spin gap (related to 
the m = plateau) survives for arbitrary 8^0 but the 
theorem is still valid because a gapless charge mode ap- 
pears simultaneously. This is the so-called Luther-Emery 
phase describing the system when H = 0. 

Plateau phase diagram - The magnetization curve for 
a fixed hole density is given in FIG. ^ Three plateaus 
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FIG. 2: Plateau phase diagram in the (H, 5) plane. The 
expected Pauli limit H p (see text) is shown in dashed lines. 
The superconducting critical field H c is also shown (the error 
bars indicate the uncertainty in determining the transition). 



are visible at m = 0,8 and 1 — S. The first one corre- 
sponds to the well known spin gap. The second one is 
expected from the YOA relation and is the most inter- 
esting case. The third one simply corresponds to the full 
saturation of the spins. In addition to these plateaus, the 
magnetization curve displays a discontinuity in its slope 
at a magnetic field H c above the m = 8 plateau. Wc 
will see that this critical field (consistent with a second 
order transition) is the superconducting critical field of 
the system. The critical fields of the transitions delim- 
iting the plateaus have been extrapolated to the ther- 
modynamic limit for various densities (see FIG. ^1 and 
are gathered in FIG. The (5 — > limit has been taken 
with 2 holes and L — > oo. The continuous behavior vs 
8 is consistent with a doping-dependent magnetization 
plateaus [ToL Hi) . The width of the m = 8 plateau is 
much smaller than the width of the m = plateau and 
vanishes for 6 ~ 0.25. In the 5^0 limit, it is exactly the 
binding energy of the hole pair-magnon bound state pre- 
viously discussed 01- This bound state originates from 
the coupling of the magnon to the charge mode which 
leads to an abrupt decrease of the spin gap in the 8^0 
limit. Therefore, the formation of the m = 8 plateau 
is based on a similar mechanism. To illustrate this, the 
local hole and spin densities in the ground state have 
been computed by DMRG and are displayed in FIG. |3| 
for increasing magnetization on a system with 6 holes. 
For the chosen parameters, holes are paired up and one 
can see three maxima in the hole density, which qualita- 
tively correspond to each pair of holes (FIG. EJ m = 0). 
When a single magnon is present in the system (FIG. 
to = 0.021), it preferably binds to a pair of holes. This 
is due to the fact that holes gain kinetic energy in a fer- 
romagnetic environment. As more magnons are added, 
more bound states progressively form until m = 8, i.e. 
when the number of magnons equals the number of hole 
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FIG. 3: (Color online) Evolution of local densities for increas- 
ing magnetization m in a system with 6 holes and L — 48. 
The spin density is normalized by 25 /m so that spin and hole 
densities always have the same mean value. 



pairs. Adding another magnon in this system costs a fi- 
nite energy gap leading to the m = S plateau. Note that 
other mechanisms could lead to plateaus in this system 
but they would correspond to different parameters in the 
Hamiltonian. By analyzing local densities beyond the 
to = S plateau, one observes that hole pairs survive up 
to a magnetization to ~ 0.22. Above the corresponding 
critical field, consistent with the value of H c in FIG. ^ 
the system behaves qualitatively like a gas of decoupled 
spinons and holons (FIG.|31 to = 0.229). 

To confirm the fact that H c is the superconducting 
critical field, we have first computed the pairing energy 
A p as a function of magnetization using 



E(n h -l,S z 
-E(n h ,S z )- 



1-1/2) 
E{n h - 



- E(n h 
2,S Z ). 



1,S Z - 1/2) 



The results plotted at fixed density on FIG. 01 show that 
pairing survives far beyond the plateau. We will see that 
the superconducting correlations are also significant up 
to H c . The H c line on FIG.|21was evaluated from points 
where A p (5 2 ) = on a system with L — 64. To estimate 
the Pauli limit in this system, we identify the condensa- 
tion energy A° at H = to A p (S z — 0). Comparing it 
with the Zeeman stabilization of two up-electrons gives 
the Pauli limit Hp = Ap 1 . Extrapolated results are shown 
in FIG.|5J The fact that the actual superconducting crit- 
ical field H c is larger than Hp suggests that the Pauli 
limit is exceeded in this system. We will see that this 
can be explained by the appearance of FFLO-like pair- 
ing above the to = plateau. 

Evolution of the ground state - In order to have a bet- 
ter characterization of the nature of the ground state un- 
der a magnetic field, we have calculated superconducting 
correlations (A^ CT , (r) (0)) where A characterizes the 
pairing channel, having a discarded weight of order 10~ 6 
with 1000 kept states. Both singlet and triplet channels 
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FIG. 4: (Color online) Pairing energy A p (see definition in 
text) as a function of H for 8 — 0.125. The dashed line 
indicates the transition to the FFLO state (q — — > q = nm). 
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FIG. 5: (Color online) (a) Superconducting correlations in 
real space showing q — Trm oscillations for rung-rung sin- 
glet and leg-leg triplet (S z = 0) correlations while rung-rung 
triplet (S z = +1) correlations decay faster (continuous lines 
are fits), (b) same in log-log scale showing the algebraic de- 
caying with exponent a (straight lines are envelopes of fits) . 



have been studied for nearest-neighbor (NN) pairs that 
can be either along a rung (called " rung" pair) or on the 
same leg (called "leg" pair). When H ^ 0, the four types 
of spin-pairing are inequivalent 



singlet: A|, (r) = 
A*(r) 



triplet: 



A*,(r) 



Cr] c r+p] and A'^(r) 



in which p determines the position of the partner (on 
same rung or same leg). In the m = plateau, two-leg 
superconducting ladders have anisotropic d-wave singlet 
pairing Rung-leg correlations and leg-leg/rung- 

rung correlations thus have an opposite sign in the sin- 
glet channel. Above the m = plateau, we found that 
both singlet and leg-leg S z = triplet channels give rise 
to dominant correlations which long-range behavior is 
cos(qr)/r a (see FIG. EJ). The wave vector of the oscilla- 
tions varies as q — irm. This can be qualitatively under- 
stood by considering the Fermi momenta of paired elec- 
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(1 — 5 ± to) on a single chain [20j. Cooper 
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FIG. 6: (Color online) Inverse of the superconducting corre- 
lations exponent a as a function of H for two densities. Sin- 
glet and S z — triplet correlations coexist above the m = 8 
plateau. Error bars are estimates of uncertainties from fitting. 



pairs with S z — are formed with a total momentum 
q = kp — kp — irm characteristic of FFLO phases Q. 
We extracted the exponent a to determine the strength 
of the correlations (see FIG.|H|). In spite of the difficulty of 
extracting precisely a, we found clear evidence that both 
channels arc coexisting above the plateau (although, be- 
cause of remaining finite size effects, we cannot conclude 
on the definite dominance of one w.r.t. the other |2cj). 
Such a rich behavior emerging under a parallel magnetic 
field might be surprising at first sight. However, a strong 
field dependence of the superconducting correlations is 
in fact expected since pairing originates from magnetic 
fluctuations. In the irrational plateau, weaker supercon- 
ducting correlations are found together with an anomaly 
in the compressibility (not shown). The latter could be 
explained from the finite energy gap needed to remove 
a pair of holes at fixed magnetization (leading to a de- 
viation from the m — 8 stability condition). Dominant 
charge density wave (CDW) correlations could then be 
expected in this plateau |2Qj. Note also that the overall 
magnitude of the triplet signal is about a hundred times 
smaller than the singlet one probably due to a more com- 
plex internal pairing structure which NN pairs only ap- 
proximate. The leg-leg triplet pairing is symmetric vs 
the exchange of chains but further investigations have to 
be done to determine the precise orbital structure. 

Experiments - Experiments in the superconducting 
phase of the ladder material SCCO are difficult because 
of high pressure and the presence of a chain subsystem. 
Transport experiments have found that the supercon- 
ducting critical field is strongly anisotropic and that the 
Pauli limit is likely to be exceeded 0, Q • We would like 
to point out that this behavior does not necessarily imply 
polarized triplet pairing at H = or spin-orbit coupling. 
Indeed, this study suggests that ladder systems should 
have a wide FFLO phase explaining the exceeding of the 
Pauli limit. Also, possible observation of irrational mag- 
netization plateaus could determine the hole doping 5 in 



ladders. Note that a NMR study proposes 5 ~ 0.1 in the 
superconducting phase [2l| . 

Conclusion - In conclusion, we show that irrational 
magnetization plateaus predicted in do appear in 
the study of the isotropic t — J model on two-leg ladders. 
Furthermore, the system has a superconducting critical 
field larger than the Pauli limit which can be explained 
by the emergence of pairing with a total finite momen- 
tum q = nm typical of a FFLO phase. This possible 
explanation for experiments does not resort to spin-orbit 
coupling nor triplet superconductivity at H = 0. In addi- 
tion, under high magnetic field, we found the emergence 
of S z = triplet pairing coexisting with singlet pairing. 
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